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The stable and metastable states of different configurations of a loop in the form of an eight is 
studied in the presence of a magnetic field. We find that for certain configurations the current is 
equal to zero for any value of the magnetic field leading to a magnetic field independent supercon- 
ducting state. The state with fixed phase circulation becomes unstable when the momentum of 
the superconducting electrons reaches a critical value. At this moment the kinetic energy of the 
superconducting condensate becomes of the same order as the potential energy of the Cooper pairs 
and it leads to an instability. Numerical analysis of the time-dependent Ginzburg-Landau equations 
shows that the absolute value of the order parameter changes gradually at the transition from a 
state with one phase circulation to another although the vorticity change occurs abruptly. 

PACS numbers: 74.60Ec, 74.20.De, 73.23.-b 
Keywords: 



Starting with the pioneering work of Little and ParkstJ 
who studied the properties of double-connected super- 
conductors (confined samples with two surfaces, for ex- 
ample rings) in applied magnetic field there has been con- 
tinued theoretical and experimental interest in those sys- 
tems. Changing the topology of the system is expected to 
have a strong influence on the superconducting proper- 
ties, in particular in the presence of an external magnetic 
field. For instance in RefsB'H the Mobuis loop was con- 
sidered and it was obtained that the phase diagram jof 
this system differs considerably from an ordinary loopcl. 

In the present paper we investigate a different topology, 
a loop having an "eight" geometry. The latter can be 
considered as a combination of two strongly interacting 
loops (see Fig.l). In comparison with two unconnected 
loops the direction and value of the current density in one 
loop of the "eight" system will depend on the direction 
and value of the current in the other loop. As a result, the 
free energy of such a system will not be merely the sum 
of the free energy of two unconnected loops. Moreover, 
we found that for certain values of the relative sizes of 
the two loops the current density in such a loop will be 
equal to zero at any value of the applied magnetic field. 

In this work we study the stable and metastable states 
of a loop in the shape of an " eight" . This geometry can be 
obtained by grapping two radial sides of a single loop and 
twisting it over 180 degrees with respect to each other. 
We assume that there is no physical contact at the cross- 
ing point of the " eight" . We can further distinguish three 
configurations for this system as shown in Fig. 2. The 
" eight" loop was modelled as a combination of two rings 
with radius i?i and i?2- We use the Ginzburg-Landau 
(GL) theory in order to study the structure of the super- 
conducting state of such a system. We neglect the finite 
width w and thickness d of the loop. This is allowed 
when the radius of the loop is much larger than w and 
if d << A, where A is the London penetration length. 
The last condition also allows us to neglect the screen- 
ing effects (self- inductance of the system). Within these 
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FIG. 1: The configuration. 



simplifications the GL equations are one-dimensional and 
using the gauge-invariant momentum p we can write it 
in the following form 

+ /(l-/'-p')=O, (la) 
3 = fp- (lb) 

where ip = /(s)e"^^'*-' is a undimensional order parameter, 
the momentum p = V0 — A is scaled in units $o/(27r^) 
(where $o is the quantum of magnetic flux) , the length 
of the loop is i = 27r(i?i -|-i?2) and the arc coordinate s is 
in units of the coherence length ^(T). In these units the 
magnetic field is scaled in the unit Hc2 and the current 
density j in jo = c<i>o/87r^A^i^. 

At first we find the dependence of the momentum of 
the "eight" loop on the applied magnetic field. From 
equation divj = it follows that the supercurrent j ~ f^p 
is constant all over the system. Our numerical analysis 
of Eqs. (la,b) shows that in the stationary state, F does 
not depend on s. As a result p = const along the loop 
and the order parameter is equal to / = ^ 1 — . 

Let us first consider the circulation of the momentum 
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FIG. 
loop. 



2: Different ID model configurations for the "eight" 



p for system [A) of Fig. 2 

(j) p - ds = 27r(i?i + B.2)p = 27m -J A ■ ds - J A • ds, 

(2) 

where the points 1,2,3,4, are at the twist point of the loop 
which in Fig. 2 are shown slightly separated from this 
point for clarity. As a result we find 
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Ri+R: 



■{n- ($1 - $2)), 



(3) 



where $1^2 = HRf 2/2 is the magnetic flux through the 
rings with radius i?i and i?2, respectively (in $0 units), 
n = § V0o?s/27r is an integer number defining the circu- 
lation (or vorticity) of the phase of the order parameter 
in the ring . From Eq. (3) follows the interesting prop- 
erty that for Ri — R2 the momentum and hence the 
current density will be equal to zero for any value of the 
applied magnetic field in system (A). Thus this system 
will have no response to an applied magnetic field and 
consequently superconductivity should be conserved up 
to very high magnetic fields. The other consequence of 
it is that there is no oscillations in the phase diagram 
H — T oi such a system. 

It is easy to show that for the systems {B, C) the de- 



pendence p{H) are, respectively, given by 
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In the limit R2 
known result 



Rl+R; 



Ri + R2 
(n- ($1 + $2))- 



(4a) 
(4b) 



we obtain from Eqs. (3,4) the well 
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which is validnfor a single double connected loop (see for 
example Ref.a). It is interesting to note that for system 
(B) the value for p is a factor {Ri/{Ri + R2)) less as 
compared to the single loop case. 

Using Eqs. (3,4) we can find the stable states of our 
systems {A — C). This is determined by the global min- 
imum of the Gibbs free energy 



G[f] 




(6) 



resulting in a minimum value of p at given H . The de- 
pendence of p{H) is a periodic function of H with period 
2/{Rj-R^), 2/Rj and 2/{Rj + Rl) for systems (A), {B), 
(C), respectively. Note that the maximum value oi Pmax 
for the thermodynamically stable state of our systems 
{A - C) is the same: 1/2 (i?i + R2). 




FIG. 3: The Gibbs free energy of the "eight" loop in the (A) 
configuration (curve 4) with parameters _Ri = 2^, i?2 = C s-nd 
the energies of single rings (with radius R = 2^ (curve 3) and 
R — (curve 1)). Curve 2 corresponds to the sum of the 
energy of these two uncoupled rings. 

The two rings in our "eight" loop are strongly cou- 
pled and as a result the Gibbs free energy of this sys- 
tem is not simply the sum of the energies of two uncou- 
pled rings plus a small interaction||term as for the case 
of two magnetically coupled ringsoQ. In Fig. 3 the de- 
pendence of the equilibrium energy of the " eight" system 
(case (A) - curve 4) and two separate rings with different 
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radii (curves 1, 3) are presented. In the same figure the 
sum of the Gibbs free energies of the two uncoupled rings 
is shown (curve 2). The dependence G{H) of the "eight" 
loop is more similar to the dependency of G{H) for the 
separate rings then with the sum of the two energies. 

Analysis shows that the state of the system for a spe- 
cific value of the phase circulation n may be stable (more 
exactly metastable) for \p\ > Pmax- This state becomes 
unstable when the second variation of the Gibbs free en- 
ergy (6) becomes equal to zero. Let's consider small de- 
viations /j <C f,Pe ^ P (where the full solution are 
f = f + fe aijdnP = p + Pe) from the stable solutions 
of our systemsQa. The Ginzburg-Landau equations lin- 
earized with respect to the small perturbations / and p 
are given by 



/,(l-3/2-p2)_2/pp, = 0, (7a) 
fp, + 2/p/, = C. (7b) 



The constant C in Eq. (7a) has to be equal to zero be- 
cause otherwise Pe ^ and /e 7^ for any / and p. Using 
Eq. (7b) we solve for p^ and substitute it in Eq. (7a) in 
order to obtain the following equation 



M6p^-2) = 0, 



which has the following solution 



fe{s) — Acos{lljs) + Bs'm{ujs), 



(8) 



(9) 



with ui = \/ — 2. Taking into account the boundary 
condition /e(0) = /e(i) we find that a state with fixed n 
becomes unstable when the following condition is fulfilled 



(10) 



Inserting i?2 = in Eq. (10) we obtain the stability 
condition of a single loop as was obtained in Ref.EI (ex- 
pressed in <I>, n and ^/Ri). In the case of large (i.e. 
min[i?i,i?2] ^ loops Pc = l/\/3 and the current den- 
sity is equal to j — pc{l—pl) = jc = 2/3\/3 - which is the 
depairing current density. Thus Eq. (10) has a simple 
physical meaning: when the kinetic energy of the super- 
conducting condensate (or Cooper pair) becomes of the 
order of the potential energy the state becomes unstable. 
The finite length of the loop starts to play an essential 
role when L is about the Cooper pair size which results 
in different values for pc- 

Using Eq. (10) it is easy to find the critical values of 
the magnetic field Hen at which the transition from state 
n to state n+1 occurs. Consider for example system (A) 
(see Fig. 2), and combine Eq. (10) with Eq. (3) we 
obtain 




Ri + R2 




(11) 



FIG. 4: Magnetic dependencies of p{H) (a) and j{H) (b) for 
the three systems shown in Fig. 2. Curves 1,2,3 correspond 
to the systems (A), (B), (C), respectively, with Ri = 2^ and 

R2=C 



It is interesting that for the configurations {B) and (C) 
the critical fields will be different but Pc (and Pmax) will 
be the same under the condition that the length of the 
loops are the same. Moreover for loops with L/S^':^ 1 the 
value of the critical momentum does not depend on the 
size of the system. Therefore, we may conclude that con- 
dition (10) is an universal condition for one-dimensional 
double-connected systems of the type shown in Fig. 2. 
Even if the loops have a different shape, for example el- 
lipsoid, the condition for vortex entry will be described 
by Eq. (10) (with the change of (i?i i?2) by L/2tt) be- 
cause the order parameter in such a sample will also be 
uniform along the loop. Indeed, a loop of arbitrary shape 
with length L is equivalent to a one-dimensional wire of 
length L with periodical boundary conditions. The pres- 
ence of a magnetic field leads to a nonuniform (uniform 
for a single circular ring) distribution of the vector po- 
tential along the wire which is compensated by the term 
in such a way that the current deiasity and the order 
parameter are uniform along the wiret^l. 

We checked Eq. (10) through a numerical solution 
of the time-dependent one-dimensional Ginzburg-Landau 



4 





(b) 



(c) 





1.5 



FIG. 5: Dependencies of the order parameter (a), the phase 
of order parameter (b), the superconducting (c) and normal 
(d) current density when switching the applied magnetic field 
from 2.b\Hc2 to 2.52Hc2- The interval (0,2) corresponds to 
the ring with radius i?i = ^ and (2, 3) to the ring with radius 
R2 = 0.5^. Curves 1, 2, 3, 4, 5, 6 are at the times 500, 
512.5, 512.95, 513, 513.5, 550 r respectively. Curves 1, 2, 3, 
correspond to times less than to which is the time at which 
the order parameter equals zero at one point on the ring (for 
our parameters to — 512. 98r). 



equations 

d_ 

ds 



(-iV-A)2^ + ^(|^|2-l), (12a) 



^-Re(V'*(-iV- A)V') 
as 



= 0. (12b) 



Here time is scaled in units of r = 47rCT„A^ (T)/c^, the 
electrostatic potential in units of c<i>o/87r^^ACT„ ((t„ 
is the noonal-state conductivity), and 7 is a relaxation 
constantli3. 

In Fig. 4 the dependency of p (Fig. 4(a)) and the 
current density (Fig. 4(b)) on the applied magnetic field 
are shovifn for the three systems (we took Ri — 2^, R2 = 
If). The magnetic field was increased linearly from H = 
to 2Hc2 during a time period of At = 2 • IO^t. The 
theoretical value of Pc for such a system is Pc — 0.593. 
From our numerical calculation we obtained Pc — 0.596. 
The difference between the numerical result and Eq. (10) 
is in the range of our numerical error. 

From our numerical solution of the time-dependent 
Ginzburg-Landau equations we find the evolution of the 
system from a state with phase circulation 2?™ to a state 
with 27r(n + 1). In Fig. 5 the dependence of \4>\, </>, js 
and jn are shown for the system {A) with parameters 
Ri — R2 = 0.5f . The magnetic field was increased 
gradually to 2.51i?c2 and then in one step to 2.52Hc2- 
At time t = the "eight" loop was at a metastable 
state at H = 2.51Hc2- At t = <o — 513t (curve 4 in 
Fig. 5) the order parameter reaches zero at a certain 
point on the ring. Simultaneously the phase difference 
A(j) = (t){+0) — Mr^) near the zero point becomes equal 
to — TT. In Refs.EJtJ the vanishing of the order parameter 
in one point on a single ring was considered and it was 
found that when such a solution exists the superconduct- 
ing current density is equal to zero and A(/) — ±tt (sign 
depends on the current direction) . Although in our time- 
dependent problem the full current is not equal to zero, 
nevertheless there is a similarity between the stationary 
and the non-stationary problems - in both cases when 
the order parameter reaches zero there is a phase shift of 
A(f) — ±7r near the point where '0 = 0. 

For t > to the superconducting current (see Fig. 5(c)) 
changes sign near the minimum point and hence it leads 
tOj-a-phase change of A(j> = 27r (from — tt to +Tr) according 
toliJ (see also Fig. 5(b)). As a result a phase circulation 
of ^ \7(j)ds=2Tr appears in the ring ( see Fig. 5(b)). Then 
the order parameter increases in the point ipmin and 
gradually decreases to (see Fig. 5(b)). Thi|S-|is-.also in 
qualitative agreement with the results of Refs.E3o where 
it was found that for increasing IV'I mm (or \js\) A(f) de- 
creased to zero near the minimum point. 

Because is a continuous function everywhere in the 
loop, the current density (superconducting js and normal 
jn) will also change continuously. The time evolution of 
js and jn are shown in Fig. 5(c) and Fig. 5(d), respec- 
tively. These currents are a function of the arc-coordinate 
during the transition process and only the sum js + jn 
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does not depend on s. For large times the normal current 
density decays and the order parameter becomes uniform 
in the system. 

In conclusion, the stable and metastable states of three 
different configurations for the " eight" loop in a magnetic 
field were studied. We found that the state with fixed 
phase circulation becomes unstable when the value of the 
momentum of the superconducting electrons reaches a 
critical Pc- Then the kinetic energy of the superconduct- 
ing condensate becomes of the same order as the potential 
energy of the Cooper pairs resulting in an instability. Nu- 
merical analysis of the time-dependent Ginzburg-Landau 



equations shows that the absolute value of the order pa- 
rameter changes gradually at the transition from a state 
with one phase circulation to another although the vor- 
ticity changes abruptly. 
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